Dilute-Bose-Gas Approach to ground state phases of 
3D quantum helimagnets under high magnetic field 

Hiroaki T. Ueda 1 , Keisuke Totsuka 1 , Tsutomu Momoi 2 

| : Yukawa Institute for Theoretical Physics, Kyoto University, Kitashirakawa Oiwake-Cho, 

! Kyoto 606-8502, Japan 

' 2 Condensed Matter Theory Laboratory, RIKEN, Wako, Saitama 351-0198, Japan 



Abstract. We study high- field phase diagram and low-energy excitations of three-dimensional 
quantum helimagnets. Slightly below the saturation field, the emergence of magnetic order may 
be mathematically viewed as Bose-Einstein condensation (BEC) of magnons. The method of 
dilute Bose gas enables an unbiased quantitative analysis of quantum effects in three-dimensional 
. helimagnets and thereby three phases are found: cone, coplanar fan and an attraction-dominant 

one. To investigate the last phase, we extend the usual BEC approach so that we can handle 2- 
magnon bound states. In the case of 2-magnon BEC, the transverse magnetization vanishes and 
long-range order occurs in the quadrupolar channel (spin-nematic phase). As an application, 
we map out the phase diagram of a 3D helimagnet which consists of frustrated J1-J2 chains 
coupled by an interchain interaction J3. 
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Magnetic frustration introduces several competing states which are energetically close to 
each other and thereby destabilizes simple ordered states. One way to compromise two 
\q ■ or more competing orders is to assume a helical (spiral) spin structure jT|. In this letter, 

we discuss the high-field behavior of a spin-1/2 Heisenberg model with generic interactions: 
H = Jij Sj-Sj + H^^- Sj. For the simplest case with one magnetic ion per unit cell, one 

can easily find the classical ground state by minimizing the Fourier transform of the exchange 
interactions: 
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where the summation is taken over all j-sites connected to the i-site by Jy. When e(q) takes its 
^ . minima e m i n at q = ±Q, helical order with the wave number Q or — Q appears (±Q are 

not equivalent to each other). When the external magnetic field is applied in a direction 
perpendicular to the spiral plane, the spiral at H=0 is smoothly deformed into the so-called 
cone state (figure Q]) as H is increased and this persists until all spins eventually get polarized 
at the saturation field H c [2]. 

One of the simplest models which exhibit, at least in the classical limit, the helical order is a 
three-dimensionally coupled Heisenberg chains with nearest-neighbor- (NN) J\ and next-nearest- 
neighbor (NNN) J2 coupling. Rise of multiferroics revives study of helimagnetism and many 
compounds which contain these ID-chains as subsystems have been reported (see, for instance, 
TABLE I. in reference [3]). For example, a helimagnetic material LiCuV04 may be viewed 
as coupled quantum S = 1/2 J\-J% chains and exhibits helical spin order, which is expected 
from the classical theories, and ferroelectricity simultaneously under moderate magnetic field. 
When the field is very high, on the other hand, this compound shows modulated collinear 



order, which contradicts with the aforementioned classical prediction, and this may suggest that 
quantum fluctuation plays an important role [4j. Therefore, it would be interesting to explore 
the possibility that quantum fluctuation replaces the classical cone state with other stable ones. 

By using dilute-Bose-gas approach, Batyev and Braginskii [5] investigated magnetic structures 
near saturation (H=H C ) and concluded that a new coplanar fan phase appears if a certain 
condition for the bosonic interactions is satisfied. Our aim in this letter is to determine the 
stable spin configurations of a specific 3D spin-1/2 helimagnet in a fully quantum-mechanical 
manner. 
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Figure 1. Two spin structures considered Figure 2. Ladder diagram of the 
here: 'cone' (left) and 'fan' (right). In the fan interaction vertex 
structure, spins are lying in a single plane. 
The Q-vector is pointing along the c-axis. 

General Formalism- By taking fully polarized state as the vacuum and treating spin-flips as 
hardcore bosons Sf = — 1/2 + pjfii, = f}J, Sf = fy, we can rewrite the spin Hamiltonian as: 

H = J>(q) - M )/3t/3 q + ^ E V*A*A?-*hfo> ( 2 ) 

1 q:k,k' 

w(q) = e(q) - e min , y, = H c - H , H c = e(0) - e min , F(q) = 2(e(q) + U) , (3) 

where U(—>oo) has been added to insure the hard-core constraint and the minimum of 
e(q) is taken e m i n = e(±Q) for helimagnets. In what follows, we consider a cubic 
lattice (we reserve (a, b, c) to label the three crystal axes) and assume that helical- and 
ferromagnetic/antiferromagnetic order occur along the c-axis and in the ab plane, respectively 
(i.e. Q = (0,0, Q) or Q = (jr, 7T, Q)). On general grounds, we may expect that magnon BEC 
occurs when the external field is H < H c (// > 0). 

The thermal potential per site E/N of the dilute Bose gas is determined by the renormalized 
interaction (t-matrix) among the condensed bosons at q = ±Q and the ground state densities 
P±Q are obtained by minimizing E/N. If we denote the renormalized interactions between the 
like bosons and that between the different ones respectively as T± and I^, the energy density 
E/N is given by 

E 1 

= 2 Tl ( p Q + P-q) + r 2 PQ P-Q ~ P(PQ + P-q)> ( 4 ) 

where /? q = |(/3q)| 2 /iV. Different phases appear according to the values of r^. When 
T2 > Ti > 0, E/N is minimal for the choice pq = p = p/Ti, p_Q = (or vice versa) and 
E/N = —p?/{2Ti). Then, the spin configuration is determined as: 

<A) = Vpexp{±i(Q.R, + #)}, (Sf) = ~ + p, 

(Sf) = Vpcos(Q-R, + 9) , {Sf) = Tx/psinCQ-R; + 9). 



That is, the cone state (the left panel of figure [I]), which exists already in the classical case [2], 
is favored for 1?2 > IV 

If Ti > T2 and r = T± + T2 > 0, on the other hand, the two modes condense simultaneously 
and the ground state is determined as: pq = p-Q = p' = p/T, -j| = —p 2 /T 

(6) 



(Sf) = -~ + 4p'cos 2 (Q-R* + HL_^) , <s±> = 2 v / 7cos(Q-R / + _J_^ )e T^ 



The existence of the two phases 9\ and 62, which correspond respectively to the two condensates 
(/3q) and (/3-q), lead to two different low-energy excitations. Since (Sf)/(Sf) = — tan ei + 6 ' 2 ; 
the spins assume a coplanar configuration (fan) shown in the right panel of figure [TJ According 
to the standard theory [6], this phase does not exhibit ferroelectricity. 

In these two Bose condensed phases, there exists two low-energy modes at q = ±Q and the 
low-energy physics is described by the effective Lagrangian with U(l)x U(l) symmetry; one 
comes from the axial (around the external field) symmetry and the other from an emergent 
translational symmetry. In the cone phase, only one of the two bosons condenses and there is 
one gapless Goldstone mode (the two U(l)s are no longer independent). Meanwhile, the fan 
phase breaks both symmetries and has two types of gapless Goldstone modes(for more detail, 
see reference [7]). 

When Ti < or Ti + T2 < 0, low-energy bosons around q = ±Q attract each other. If 
the energy ^) is taken literally, one may expect, on general grounds, first order transitions 
to occur (provided that the cubic terms of p are positive). However, eq.([2]) is based on the 
assumption that magnon BEC occurs in the single-magnon channel and may not work when we 
expect magnon bound states stabilized by strong attraction. In fact, according to the standard 
Bethe-Salpeter method, divergingly large (renormalized) attraction implies the existence of 
stable 2-magnon bound states as a pole in the two-particle Green's function corresponds to a 
bound state. Therefore, the simple single-magnon BECs will give way to those in multi-particle 
channel. To check the above scenario, we solve the 2-magnon scattering problem. We denote 
the ladder diagram as T(A,K;p, p') (see figure [2]), where we take the total energy of initial 
state as —2p — A. Then, on the fully saturated ground state, the exact scattering amplitude 
M(A, K; p, p') = r(A, K; p, p') + T(A, K; p, — p') is obtained by solving the following integral 
equation: 

M(A,K,p, P ) = V(p - P )+n-P -P)-^vL w (k/2 + p")+ W (k/2-p") + A-iO+ " 

(7) 

Even in the presence of the boson condensate, as long as it is dilute, we can safely use 
this scattering amplitude. Hence, one obtains the renormalized interaction Ti and T2 as 
ri = (1/2)M(0,2Q;0,0), T 2 = M(0,0;Q,Q) 0. Moreover, the pole of M (A, K; p, p') gives 
the stable bound state. 

Coupled J1-J2 model- Having established the formalism, we proceed to investigating a specific 
model-a frustrated spin- 1/2 model on a simple cubic lattice whose Hamiltonian is given by 

H= {«JlSr-S r +g + J2S r -S r +26 c + JsSr'Sr+gj} + H^^SJ, (8) 

r,i=a,b j 



In (|8|), the J1-J2 chains are running in the c-direction and J3 controls the coupling among 
adjacent chains. The wave number Q characterizing the condensate is given either by Q = 
(0,0,(5) (^3 < 0) or by Q = (tt, 7r, Q) (J3 > 0) where Q = arccos(— J1/4J2). Hence the spiral 



occurs in the c-direction. To determine the spin structure of our J1-J2-J3 model, we solved 
eq.© by assuming the following form [5]: M(A, K;p, p') = (M) + Aicosp' c + A2Cos2p' c + 
Ascosp'a + A^cosp' b , where (M) and A{ are functions of A,K and p. After some numerical 
calculations, we obtained the phase diagram shown in figure El We show only the frustrated 
region —4 < J1/J2 < 4, where cone structure with incommensurate Q is expected classically. 

Figure 3. Phase diagram slightly below 
saturation (H < H c ) mapped out in (Ji, J3)- 
plane (J2(> 0) is used to set the energy 
unit), (i) cone phase and (ii) coplanar fan 
phase are one-magnon condensed phases. In 
the phase- (iii), attractive interaction may 
imply instabilities toward other phases e.g. 
conventional ferromagnetic one or more exotic 
multipolar ones. The phase (iv) (nematic) is 
characterized by the condensation of the 2- 
magnon bound state and leads to the nematic 
order in the transverse direction. The region 
I J3 \/J\ 1 is omitted for a reason described 
in the text. Inset: The same phase diagram 
for the large negative interchain coupling 
(-J3 > 0.1). 

For J3 — > 0, low-energy quantum fluctuation destabilizes T and our approach cannot be 
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extended to J3 = continuously (T becomes 0(J S ) and Ti — > T2 at the leading order in J3). 

Finally, we compare the phase diagram shown in figure [3] with that of the ID J1-J2 chain 
(J3 = 0). Let us begin with the case J\ > 0. Near saturation, two dominant phases are found in 
ID [9j: (i) 'chiral phase (VC)' with finite vector chirality parallel to the magnetic field and (ii) 
'TL2' phase where the system is described by two Tomonaga-Luttinger (TL) liquids. Obviously, 
the former turns, after switching on J3, into the cone phase. A close inspection of the two 
gapless TL modes near saturation tells us that the TL2 phase should be identified with the fan 
phase here. 

For the ferromagnetic case J\ < 0, BECs of n-bound magnon states (n > 2) are expected in 
ID chain jlOj . We found how strong inter-chain coupling destabilizes the 2-bound magnon BEC. 
The study of the stability of higher n-bound magnon state (n > 3) in 3D is a future problem. 
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